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ABSTRACT 

In  a recent  paper,  Proschan  and  Singpurwalla  (1978)  have  proposed  a 
new  approach  for  making  inferences  from  accelerated  life  testa.  Their 
approach  is  significantly  different  from  those  that  have  been  considered 
in  the  past , and  is  motivated  by  what  is  actually  done  in  practice . A 
priori  information  which  is  generally  available  to  the  engineer  is  incor- 
porated under  their  approach  by  adopting  a Bayesian  point  of  view.  The 
usual  assumptions  about  the  failure  distributions  and  the  acceleration 
functions,  which  are  appealing  from  a statistical  point  of  view,  have 
been  sacrificed  to  achieve  greater  reality. 

In  this  paper,  we  apply  this  new  approach  to  some  real  life  data 
arising  from  an  accelerated  life  test.  In  the  process,  we  explain  our 
new  approach  in  a manner  which  makes  it  easy  to  understand  and  apply  by 
the  reliability  practitioner. 
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1*  Introduction . 

In  practice,  long  life  items  are  subjected  to  larger  than  normal 
stresses  in  order  to  obtain  failure  data  in  a short  amount  of  test  time. 
Such  tests  are  called  accelerated  or  overstress  life  tests , and  the  goal 
is  to  make  inferences  about  the  life  distribution  ‘of  the  items  at  the 
normal  stress  levels  using  failure  data  from  accelerated  tests. 

The  current  approach  to  this  problem  involves  making  assumptions 
about  the  distribution  of  failure  times,  and  about  the  functional  rela- 
tionship betveen  the  parameters  of  the  failure  distribution  and  the  applied 
stress.  Such  a relationship  is  known  as  an  acceleration  function  or  a 
time  transformation  function ; examples  of  these  are , the  Power  Law , the 
Arrhenius  Law,  the  Eyring  Law,  etc.  . Another  assumption  commonly  made 
states  that  at  all  stress  levels , the  failure  times  are  governed  by  dis- 
tributions which  sure  members  of  the  same  parametric  faimily,  such  as  an 
exponential  or  a Welbull.  These  assumptions,  though  appealing  from  a 
statistical  point  of  view,  may  be  unreasonable  in  many  practical  situations. 
Of  particular  concern  is  the  last  assumption.  The  different  stress  levels 
may  have  different  effects  on  the  "failure  mechanisms" , and  thus  from  an 
engineering  point  of  view,  it  may  be  more  realistic  to  allow  for  differ- 
ent forms  of  the  failure  distribution  at  the  different  stress  levels. 

In  a recent  paper,  Proschan  and  Singpurwalla  (1978)  (hereafter  re- 
fereed as  PS(1978))  have  proposed  a new  approach  for  making  inferences 
from  accelerated  life  tests  which  requires  neither  distributional  assump- 
tions nor  the  specification  of  a time  transformation  function.  Rather, 
their  approach  is  Bayesian,  and  is  motivated  by  procedures  actually  used 
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in  practice  by  reliability  engineers.  The  Bayesian  point  of  view  enables 
a user  to  incorporate  sene  a priori  information  which  is  available  in 
accelerated  life  tests. 

In  Section  2 we  shall  review  the  pragmatic  Bayesian  approach  to  ac- 
celerated life  testing  by  PS(ip78),  and  present  it  in  a manner  which 
makes  it  easy  to  understand.  In  Section  3 we  shall  demonstrate  the  use- 
fulness of  this  approach  by  applying  it  to  some  real  life  data  on  accel- 
erated life  tests  presented  in  Nelson  (1970). 

Our  objective  in  writing  this  paper  is  to  make  this  new  approach  to 
accelerated  life  testing  accessible  to  the  reliability  engineer,  and  to 
demonstrate  its  usefulness  by  applying  it  to  a realistic  situation. 
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2*  A Pragmatic  Bayesian  Approach  to  Accelerated  Life  Testing. 

We  shall  denote  the  k accelerated  stresses  (environments)  by  ' 

• • • « ^ the  normal  or  use  conditions  stress  (environment) 
by  Eu.  Let 

E1  > E2  > ...  > ^ > Eu,  (2.1) 


where  > E^  denotes  the  fact  that  E^  is  more  severe  than  E^ . 

Let  Fj  be  the  failure  distribution  of  the  items  tested  under  environ- 
ment Ej , and  let  A^(x)  denote  the  failure  rate  of  Fj  at  time  point  x 2 0. 

Because  of  condition  (2.1)  it  is  logical  from  a physical  point  of 
view  to  assume  that  for  any  x a 0, 


A^(x)  a Ag(x)  a ...  a A^(x)  ^ A^(x)  • (2.2) 

Using  failure  data  obtained  under  E^,  Eg,  ... , E^ , we  would  like  to  obtain 
*j(x),  an  estimate  of  Aj(x),  J « 1,  2,  ....  k,  such  that  for  Booe 
0 i L < •,  and  all  x e [0,  l], 

* . rt  * t st  st  » 

A^x)  a Ag(x)  a ...  a A^x).  (2.3) 

The  notation  X ^ Y denotes  the  fact  that  X is  sto<>>»<mtiealiy 
than  T;  that  is 

P[X  a x]  a P(Y  a x]  for  all  x. 

In  order  to  obtain  estimators  Aj(x) , J ■ 1,  2,  ...,  k,  which  satisfy 
(2.3) » we  shall  use  a Bayesian  approach. 


Under  this  approach,  condition 


2 


Let  Nj(t)  be  the  number  of  items  tindergoing  life  test  at  time  t 

under  environment  E^.  Divide  the  time  interval  [0,  L]  into  sub-intervals 

of  length  h > 0,  where  h is  chosen  to  make  L a multiple  of  h.  Thus , the 
total  number  of  sub-intervals  is  L/h. 

For  convenience,  we  shall  Introduce  the  following  notation: 


tA  ■ the  time  point  (i  - l)h,  i ■ 1,  2,  (L/h). 


[t^,  t^  ♦ h)  ■ the  i^htime  interval 


Let 


B,  , ■ the  number  of  items  exposed  to  the  environment  E at  time  t, 

J t1  J 1 


X,  . ■ the  failure  rate  of  F,  at  t. 
j J i 


,th 


x.  ^ ■ the  number  of  failures  in  the  i time  interval  under 
J ’ environment  Ej 

p.  . ■ the  probability  of  failure  of  a unit  in  the  ith  time  interval 
J ’ under  environment  Ej. 


PML 


-i-i 

lPj 


then  Is  the  conditional  probability  that  an  item  which  has  survived 
to  time  t.  will  fall  by  time  t.  will  fall  by  time  t.  + h.  We  can  also 
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interpret  p*  ^ as  the  average  failure  rate  over  the  interval  [t^,  t^  + h). 

If  we  assume  that  there  are  no  withdrawals,  removals,  or  censoring,  then 


NJ,i+l  * NJ,i  " xJ,i* 


i - 1,  2 (L/h). 


Our  Bayesian  analysis  involves  assigning  prior  distributions  to  the 

average  failure  rates.  Suppose  that  the  p*  , i * 1,  2,  ...,  (L/h),  are 

J f* 

a priori  independent  beta  random  variables  with  prior  parameters  a and  $ , 

J 

ao  that  their  marginal  densities  are 

r(o  + n i &.-1 

f(pJ,i)  * r(o)r(Bj)  (pJ,i)  (1  “ pJ,i) 

Then,  given  the  N,  . 's  and  the  x,  . 's,  the  posterior  density  of 
J **  J 9* 

(P3.i’PJ,2’  -•Pj;L/h))  18  [cf*  PS(1978)1 


(L/h) 


r(a  ♦ 6,  ♦ H,  ,) 


e.+N,  ,-x.  4-l 


• Mj  ' s j/  . T j 4-* 


2 .2 . A Bayesian  Estimation  of  the  Ordered  Average  Failure  Rates. 

In  the  Bayesian  context,  condition  (2.2)  leads  us  to  the  requirement 
that  for  every  fixed  value  of  i , p*  . 1*  p*  . for  J » 2 , 3 , 

J— X »1  J >1 

Thus , the  prior  distributions  of  p*  . and  p*  will  have  to  be 

J—A  ,1  J ,1 

chosen  so  as  to  ensure  that 

PlpJ_ltl  * p]  i plP*fl  * p]  for  all  p,  0 s p s 1.  (2.5) 
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One  way  of  achieving  the  above  condition  is  to  require  [see  PS(1978), 
Appendix  A]  that 


6j  * for  J = 2,  3,  . ...  k. 


To  be  assured  that  Condition  (2.5)  is  also  satisfied  with  respect 

to  the  posterior  distributions  of  p*  , J and  p*  , , it  is  sufficient  that 

J— J-  >1  J ,1 

for  every  fixed  value  of  i. 


x i x 

J,i  *J-l,i 


and 


(2.6) 


ej+NJ,i-xJ,i2ej-l+NJ-14"Xj-l,i 


for  J = 2,  3,  ....  k. 

■^e  first  part  of  condition  (2.6)  states  that  the  number  of  failures 
in  the  interval  [t±,  t±  + h)  under  environment  E.  must  not  be  greater  than 

the  number  of  failures  in  the  same  interval  under  environment  E , for 

J 

all  values  of  J.  Furthermore,  as  is  discussed  in  PS(l978),  a reasonable 
strategy  is  to  have  more  items  initially  on  test  under  the  more  severe 


environment  E,  ^ than  under  the  environment  E. , so  that  N,  < N 
for  all  J. 

Since  ■ Hj the  second  part  of  condition  (2.6) 
can  be  written  as 


0J-1  + NJ-l,i+l  S 


B,  + H 


J,i+1 


for  J - 2,  3, 
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Consider  the  time  interval  [(i-l)h,  ih) ; by  Condition  (2.6)  we  require 


Xl,i  * X2,i  * * * XJ-l,i  * XJ,i  * •"  4 \,i 


+ Nl,i  " x!  i 5 e2  + N2  i " x2  i ^ •••  5 + Nj  i i " xj  i j 


jt  ’ ”k,i  - \i- 


If  a reversal  occurs,  that  is,  if  either 


XJ-l.i  < XJ,i 


or  if 


6J-1  + " XJ-l,i  > + NJ,i  _ XJ  ,i  * 


then  we  pool  the  violators  and  replace  them  as  shown  below. 

Replace  both  x.  , , and  x,  . by  i(x,  . , + x,  .)  and 

“,<i  8j-i  * “j-i.i  - xj-i,i  “a  * "j.i  - xj  ,i w 

a<BJ-l  * V Vl.i  * "j,l  * XJ-1,1  ‘ XJ,i)' 

'/e  now  test  the  new  sequences  to  see  if  they  are  properly  ordered'. 
If  not,  we  replace  the  adjacent  violators  by  their  appropriate  averages. 


Thus  if 


2*xJ-l,i  + XJ,i)  " 2^xJ-l,i  + xJ,i)  < xJ+l,i* 


then  we  replace  each  of  the  three  by  the  average 


3*XJ-l.i  + XJ,i  + xJ+l,i)# 


The  ssae  procedure  is  vised  for  the  (0  ♦ N , - x,  . )'s. 

J J J 


We  follow  the  above  procedure  for  all  time  intervals  and  continue 
until  all  reversals  are  eliminated.  Note  that  excessive  pooling  will 
occur  if  the  relationship  (2.2)  does  not  actually  hold,  or  if  the  environ- 
mental conditions  are  too  similar  to  each  other. 

Assuming  a squared  error  loss,  before  pooling,  the  Bayes  estimator 
pJ,i  ls  tsee  "3(1978)) 


a + x 


'J.i 


jLiL 


a + 6 . + N . , 
J J ,i 


(2.7) 


If  we  have  done  6ome  pooling,  then  the  x,  . 's  and  the  (0  + N.  . )'a  are 

J J J »i 

replaced  by  their  appropriate  pooled  averages. 


2.U.  A Model  for  Extrapolation  to  Use  t Stress. 

V.  shell  uss  the  ;*ys  to  estimate  pj_2 P^.a/h) . ‘he 


average  failure  rates  under  the  use  conditions  environment  E . 

u 

In  the  absence  of  any  physical  or  engineering  knowledge  about  the 
relationship  between  the  average  failure  rate  and  the  corresponding  stress, 
we  shall  postulate  the  following  simple  but  reasonable  relationship 
among  the  average  failure  rates. 

For  each  value  of  i,  i **  1,  2,  ...,  (L/h) , 


+ Vk-i,i + •••  - wk-ip!  ,i » 


(2.8) 


"0 

where  wQ , w1 , . . . , wk_^  are  unknown  constants . 

The  above  relationship  states  that  the  average  failure  rate  over  a 


articular  time  interval  under  the  environment  ls  a weighted  sum  of 

the  average  failure  rates  over  the  same  time  interval  under  the  environ- 

\_2*  ....  Er 
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* A 

Let  w0’  wl»  the  le«at-8quarep  estimators  of  vQ , w,  , ...»  w, 

these  can  be  obtained  routinely  from  Equation  (2.8),  but  with  > in  place 
of  p*  throughout.  Thus  for  i - 1,  2,  (L/h),  we  have 


Pu,i  “w0  +Vk,i  + •**  + vA#1 


(2.9) 


as  the  estimators  of  the  average  failure  rate  under  E . As  a consequence, 

» (L/h)*. 

we  also  have  p*  ,T  , ■ 1 _ Y «• 

^u,(L/h)+l  A £ Pu.1* 

An  estimator  of  P^t),  the  probability  of  an  item  surviving  to  time 
t under  Eu,  the  use  conditions  stress,  is 

» t/h  « 

W * i'1  - 

A 

where  the  p* ^ are  given  by  (2.9). 
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b.  An  Illustrative  Real  Life  Example. 

S 

In  this  section  ve  shall  apply  the  methodology  discussed  in  the  previous 
section  to  some  accelerated  life  test  data  given  by  Nelson  (1970).  These 
data  represent  t!  times  to  breakdown  of  an  insulating  fluid  subjected 
to  elevated  voltage  stress  levels.  For  convenience,  ve  shall  consider 
here  only  U accelerated  voltage  levels  - 36,  3U,  32,  and  30  kilovolts 
(KV's).  The  failure  times  (in  minutes)  under  the  various  stress  levels 
ere  given  in  Table  1.  Nelson's  original  data  correspond  to  7 different 
stress  levels,  but  some  of  these  contain  very  few  failure  times  and  are 
therefore  omitted. 

We  shall  assume  that  the  use  conditions  stress  is  28  KV,  and  apply 

\ 

our  approach  to  estimate  the  failure  distribution  of  the  breakdown  times 

’ 

at  this  stress. 

Following  the  notation  of  Section  2,  we  shall  choose  L to  be  100  min- 
utes, malting  the  total  number  of  time  intervals  equal  to  100.  For  our 

prior  parameters  ve  shall,  following  the  discussion  in  Sections  (2.1)  and 

■ 

^2.2),  choose  a * 1,  with  ^ ■ 1,  Bg  ■ 2,  6^  * 7,  and  3^  * 12.  The  above 

choice  of  our  prior  parameters  3^  Bg,  and  is  motivated  by  an  in- 
spection of  the  failure  times  in  Table  1.  Observe  that  even  though  the 
stress  levels  decrease  successively  by  only  2 KV's,  the  corresponding 
change  in  the  failure  times  from  one  stress  level  to  the  next  lower  ap- 
pears to  be  more  drastic.  This  may  be  due  to  a change  in  the  failure 
■echanism  as  we  go  to  the  lower  stresses.  Note  that  N,  - 15,  N„  ■ 19. 

»3  - 15,  and  - 11. 

After  computing  the  x.  's  and  the  (fi.  ♦ N . . - x,  . ) 's  for 

» •*  J J **  J ,i 

i 


12 


J“1»2,3»U,  1*1,  2,  ....  100,  and  after  pooling  the  adjacent  rlo- 

A 

lators  as  discussed  In  Section  2.3,  we  obtain  the  Bayes  estimators  p* 

J 


using  Equation  (2.7).  These  values  are  given  in  Table  2,  wherein,  for 
convenience,  we  have  limited  ourselves  to  time  intervals  of  width  5.  Ob- 
serve that  the  pooling  scheme  ensures  that  for  each  time  interval  the 


the  p*  j'  4 decrease  with  increasing  stress  level.  For  example , P*  ,20  ■ *333:., 

P*,20  * *1250»  and  *3,20  “ 811(1  Pj  " 


» , • \ 
‘ J* 


Our  next  step  is  to  use  the  entries  in  Table  2 to  obtain  estimates  of 
the  w's  given  in  Equation  (2.8).  If  we  designate  the  30  KV  level  as  the 
level  k(*U)  of  Equation  (2.8),  and  the  32,  31*,  and  the  36  KV  levels  cub 
k - 1,  k - 2,  and  k - 3,  respectively,  then,  the  lest  squares  estimation 
of  the  w's  gives  us  for  i * 1,  2,  ....  100, 


pj^  - .2360  + .18861  p*^  + .13538  p*^  + .09857  pj^. 

We  shall  now  use  the  following  equation  to  estimate  p#  . , i * 1,  2,  ...,100, 

u,i 

where  u(*5)  denotes  the  use  condition  stress  level  of  28  KV: 


p*^  - .2360  + .1886  pj^  + .13538  p*^  + .09857  p*ti. 


An  estimate  of  the  survival  probability  at  use  conditions  voltage  of 
28  KV  is  given  by 

F(t)  - 5 (1  - p*  ,). 
i-1  5,1 

A 

The  values  of  Fu(t),  t ■ 5,  10,  100,  are  given  in  Table  3,  and 


plotted  in  Figure  1 


Table  1 


i 

* 


Tlmea  to  Breakdown  of  an 
Insulating  Fluid  (la  Minutes ) Under 
Various  Values  of  the  Stress  (In  Kilovolts). 


Stress 


36  KV 

3U  KV 

32  KV 

JO  KV 

.35 

.19 

• 27 

7.7U 

.59 

.78 

.1*0 

17.05 

.96 

.96 

.69 

20.U6 

.99 

1.31 

.79 

21.02 

1.69 

2.78 

2.75 

22.66 

1.97 

o.lS 

3.91 

1*3. 1*0 

2.07 

U.15 

9.88 

1*7.30 

2.58 

U.67 

13.95 

139-07 

2.71 

U.85 

15.93 

11*1*.  12 

2.90 

6.50 

27.80 

175.88 

3.67 

7.35 

53.21* 

19l*.90 

3.99 

8.01 

82.85 

5.35 

8.27 

89.29 

13.77 

12.06 

100.58 

25.50 

31.75 

215.10 
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Table  2. 


Estimated  Values  p*  . of  xhe  Average 

— 


Failure  Rates  Under  Various 
Values  of  the  Stress  (In  Kilovolts) 


Stress 


Time  Interval 


3 6 KV 


3h  KV 
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ALSTRACT 

In  a recent  paper , Pros chan  and  Singpurvalla  (1976)  have  proposed  a new  ap- 
proach for  making  Inferences  from  accelerated  life  tests.  Their  approach  is 
significantly  different  frost  those  that  have  been  considered  in  the  past,  and 
in  motivated  by  what  is  actually  done  in  practice.  A priori  Information  which 
is  generally  available  to  the  engineer  is  incorporated  under  their  approach  by 
adopting  a Bayesian  point  of  view.  The  usual  assumptions  about  the  failure 
distributions  and  the  acceleration  Auctions,  which  are  appealing  from  a sta- 
tistical point  of  view,  have  been  sacrificed  to  achieve  greater  reality. 

In  this  paper,  we  apply  this  new  approach  to  some  real  life  data  arising 
from  an  accelerated  life  test.  In  the  process,  v*  explain  our  new  approach  in 
a mano  r which  makes  it  easy  to  understand  and  apply  by  the  reliability  practitioner 


